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RESEARCH MEMORANDUM 


AN EXPERIMENTAL INVESTIGATION OF THE UNSTEADY LIFT 
INDUCED ON A WING IN THE DOWNWASH FIELD OF 
AN OSCILLATING CANARD CONTROL SURFACE 
By David E, Reese, Jr. 

SUMMARY 


The results of an experimental investigation of the unsteady lift 
induced on a wing in the downwash field of an oscillating canard control 
surface are presented. The in-phase and out-of-phase components of the 
lift and their respective centers of pressure were measured at frequencies 
from 10 to 110 cycles per second over a Mach mmiber range of 0.6 to 0.8 
and 1.4 to 1.9 at angles of attack of 0°, 5°^ and 10°. 

The results indicated that existing theories provide a reliable guide 
for the estimation of the magnitude of the lift derivatives and centers 
of pressure at low values of reduced frequency and low angles of attack. 

An estimation of the effect of frequency on the lift derivatives was 
developed using a simple indicial function for the wing based on the lag- 
in-downwash concept. The trends of the data with frequency predicted by 
the theory were, in general, borne out by experiment at low angles of 
attack. 


INTRODUCTION 


With the present trend in airframe configurations for guided missiles 
toward low-aspect-ratio wings and tails, the contribution of wing-tail 
interference effects to the over-all stability of the configuration has 
become of major importance. An accurate appraisal of these interference 
forces is necessary not only for the usual stability calculations but also 
in the design of automatic control equipment. 

In an analysis of the response of a missile to various guidance com- 
mands, both static and dynamic forces and moments acting on the missile 
must be considered. The requirement of rapid response to these commands 
has focused increased attention on the dynamic quantities involved in the 
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calculations. One particular problem for canard-type missiles that has 
brought to light the need for information on dynamic interference forces 
is associated with the longitudinal dynamic stability of these missiles. 

In order to provide sufficient damping of the longitudinal motion of the 
missile^ damping is often supplied artificially by actuating the forward 
control surfaces in proportion to the pitching velocity of the airframe. 
Since it is possible for a large portion of the total moment created in 
this manner to arise from the interference lift on the rear surfaces^ it 
is necessary to know these dynamic interference effects with reasonable 
accuracy in order to predict the response of the missile. A knowledge 
of not only the magnitude of the interference lift but also its dependency 
on frequency is required. 

Some theoretical estimates of static wing-tail interference forces 
have been made by Nielson^ Kaattari, and Anastasio (ref. l) and Schindel 
and Durgin (ref. 2). An experimental investigation of these forces and 
moments for several wing-tail combinations is presented in reference 3 
by Schindel and Durgin. Some theoretical values for the time -dependent 
forces are given by Miles (ref. 4 ) and Martin, Dieterich, and Bobbitt 
(ref. 5) i*or very low frequencies. However, no published data on the 
experimental determination of these unsteady interference forces and 
moments are known to be available. The present experiment was, therefore, 
planned to meet this need. 


SYMBOIB 




■^ip 

kip 

M 


lift derivative in phase with control-surface position, positive 

57 • 3^jp 

upward, , per radian 

IpV^S&o 

lift derivative out of phase with control surface position (in 
phase with control surface velocity) , positive upward 

^7 3kp ^ pgj. j-adian 

lift in phase with control position, lb 

lift out of phase with control position, lb 

V 

Mach mmiber, r — r, dimensionless 

^ speed of sound ^ 

total area of wing, sq ft 
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V 

c 

f 

k 



a 

5 

6 

So 

So 

P 

9 

CO 


free-stream speed, ft/sec 

wing chord, ft 

frequency of oscillation 

reduced frequency, dimensionless 

reduced frequency at wind-tunnel resonance, dimensionless 
distance from control surface centroid to wing centroid, ft 
dynamic pressure, Ib/sq ft 

time , sec 

center of pressure of in-phase lift, positive behind wing 
leading edge, dimensionless 

center of pressure of out-of-phase lift, positive behind wing 
leading edge, dimensionless 

angle of attack of wing and body, deg 

instantaneous control-surface deflection, deg 

instantaneous control-surface angular velocity, deg/sec 

maximum control-surface deflection, deg 

maximum control- surface angular velocity, w5q, deg/sec 

mass density of air, slugs/cu ft 

phase angle between lift and angular position of control surface, 
deg 

angular frequency, 2atf, radians/sec 

THEORY 


The purpose of this investigation is to determine the unsteady lift 
induced on a wing in the downwash field of an oscillating canard control 
surface. In this section, methods will be presented from which results 
can be obtained for the theoretical values of this lift. Frequency 
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effects on the phase relationship between the induced lift and the con- 
trol deflection are of primary interest. For purposes of computation, 
the induced lift at frequencies near zero will be treated first, however, 
and then the effects of frequency will be studied. As is usually done in 
dynamic stability calculations, the amplitude of the unsteady lift will 
be divided into two components, one in phase with the control surface 
position and the other in phase with control surface velocity, the latter 
being commonly known as the out-of -phase lift. The model used for the 
theoretical analysis consisted of a canard control surface and a wing 
alone, that is, a bodyless model. The degree to which the body effects 
are considered in this report will be mentioned in a later section. 


Frequencies Near Zero 


The amplitude of the in-phase lift derivative and its center 

of pressure at values of k approaching zero were calculated by the method 
presented in reference 1. Since this procedure is commonly used in the 
calculation of the static lift and center of pressure of wing-body-tail 
combinations, it need not be discussed further here. 

The calculations for the amplitude and center of pressure of the out- 
of -phase lift derivative were based on the lag-in-downwash concept 

appearing in reference 6. In this approximation the assumption is made 
that the downwash at the rear surface will be that created by the forward 
surface at a small time z/V earlier, where Z is taken to be the dis- 
tance from the centroid of area of the control surface to the centroid of 
area of the wing. Thus 


Agw 


d6 dt 


At 


where 

Agv/ change in downwash at the wing 
d€e change in downwash at the control 
then 


Agv 


dec . I 
" d6 ^ V 


( 1 ) 


( 2 ) 


The lift on the wing due to this change in downwash can he written 


CONFIDENTIAL 


NACA RM A55F01 


C01fi5)EWTIAI* 


5 


lOP = - q^ew(Cl^)^ 

= ,ss I (3) 

Where (Ci^)^ is the lift-curve slope of the wing alone in uniform flow. 
But dec(CLa)^ is the negative lift on the wing due to a small change in 
downwash dec* This negative lift divided hy d& is then - so 

lOP = - qS6 I Ci^ (4) 

In coefficient form this lift is 

Ci^ = - 2 ~ Ci^ (5) 

which corresponds to the answer obtained hy Martin, Dieterich, and Bobbitt 
in reference 5 - 


Frequency Effects 


The effects of frequency on the lift derivatives at a Mach number of 
1.7 were calculated using the response of the wing to an impulsive deflec- 
tion of the control surface. This indicial response of the wing was cal- 
culated by methods presented by Tobak in reference 7 . The indicial 
response was then converted to frequency response through the use of the 
Duhamel Integral. 


Since the calculations necessary to define the indicial lift response 
for the present problem were quite lengthy, they were carried out only for 
a Mach number of 1.7 angle of attack of 0^. However, it was found 

for this case that a good approximation to the frequency response obtained 
by the more exact solution could be made by assuming an indicial lift 
response based on the lag-in-downwash concept. This is an extension then, 
to include the effects of frequency, 
of the theory on which the low- 
frequency values of the out-of-phase 
lift derivative were based as given 
in the preceding section. In this 
approximation it is assumed that the 
wing experiences no lift until the 
downwash Impulse from the control 
surface has reached the midchord of 
the wing. At that time the lift of 
the wing immediately assumes its 
steady-state value. This results in 
an indicial function of the form 
shown in sketch (a), where t* 



IS 


Sketch (a) 
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the time required for the pulse to travel from control surface midchord 
to wing midchord and Cl* is the steady-state value of the lift of the 
wing in the presence of the downwash from the control. 


The time t* has an important physical meaning as is shown in refer- 
ence 7. In that paper it is shown that the indicial lift of the wing due 



Sketch (h) 

plified indicial response and those 
indicated. 


to the downwash from the suddenly 
deflected control surface has the form 
shown in sketch (h) . It is also shown 
that, with the fulfillment of certain 
mild conditions, the time at which the 
indicial curve crosses the zero lift 
axis is veiy nearly equal to the time 
t* defined above. It was also found 
that, in frequency-response calcula- 
tions of the derivatives Ct^ and 

-*-'0 

the effects of the two shaded 
peaks tended, in large part, to 
cancel. Thus good correspondence 
between the results based on the sim- 
based on the more exact solution is 


Once the indicial function has been defined, it is possible to obtain 
the variation of the lift on the wing with frequency through the use of 
the Duhamel integral. In reference 7 Tobak has used the Duhamel integral 
to derive the following equations for the in-phase and out-of-phase lift 
derivatives in teniis of the indicial functions of the wing: 


( 6 ) 
( 7 ) 

where 

Fi(e) = ci^(e) - 
FsCe) = Cl^( 0) - CLq* 


CLg = Fi(0)sin k’0 d0 + j F2(0)cos k’0 d.0 

poo poo 

= CLq^* - / Fi(0)cos k*0 d0 - kM F2(0)sin k«0 d0 


and 

CLg(0) indicial response of the wing to a step deflection & of the 
control surface 
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steady-state lift on the wing due to a deflection 6 of the 
control surface 


ciq(e) 


indicial response of the wing to a step change in pitching velo- 
city q of the control surface 


e 

k» 


steady- state lift on the wing due to a pitching velocity q of 
the control surface 

nondimensional time parameter = 

reduced frequency, defined in reference 7 


It can he shown that, for the problem considered in this section, 
and the integrals containing F2(0) are small relative to the 
remaining terms in the equations. Making this approximation and noting 
that k* = y k, we reduce equations (6) and (t) to 


C 






sin Y k0 do 


(8) 


X OO 

F3_(0)cos ^ k0 d0 


(9) 


Performing the integrations indicated in equations (8) and (9) on the 
indicial function defined in sketch (a), we derive the following relation- 
ships ; 


cos 2k 



I „ ^ sin 2k 

c ~k~ 


( 10 ) 

( 11 ) 


These expressions give the values of the in-phase and out-of-phase lift 
derivatives as functions of frequency. It is apparent that the limiting 
values of equations (lO) and (ll) as k approaches zero give the results 
obtained in the preceding section. 

The phase angle between the induced lift and the control deflection 
is 


cp = tan 


“1 = 
Lip 


tan 


-1 


(coc/2V)Ci^ 




( 12 ) 
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Substituting the expressions from (lO) and (ll) in equation (12), we find 
that 


In order to obtain the variation of the centers of pressure of the 
in-phase and out-of-phase lifts with frequency it was necessary to obtain 
values for the in-phase and out-of-phase pitching moments corresponding 
to those found for the lift in equations (lO) and (ll). For this problem 
the indicial variation of the pitching moment with time was again that 
shown in sketch (a). While the calculations are not shown, it is apparent 
that the expressions for the pitching moments will be identical to equa- 
tions (lO) and (ll) with the substitution of the appropriate moment deriva- 
tives for the lift derivatives. Since center of pressure is defined as 
pitching moment divided by lift, it can be seen that both the centers of 
pressure of the in-phase and out-of-phase lifts are equal to the center of 
pressure for the steady-state case and that there is no variation of center 
of pressure with frequency. Therefore the method of reference 1 was used 
to calculate these values. 


The investigation was conducted in the Ames 6- by 6-foot supersonic 
wind tunnel which is of the closed-return, variable-density type and which 
has a Mach number range of 0.6 to 0.9 and 1.2 to 1.9. A detailed descrip- 
tion of the flow characteristics of this wind tunnel can be found in 
reference 8. 


The model consisted of a rectangular wing and a rectangular control 
surface mounted in canard arrangement on a slender cylindrical body of 
fineness ratio ll<-. The dimensions of the model are given in figure 1. 

Both the wing and the control surface were composed of 5-percent-biconvex, 
circular-arc airfoil sections. A 5-inch, hollow, circular cylinder fitted 
with an ogival nose served as the body. The model was fabricated from 
steel with the exception of the nose and wing. The nose section was made 
of fiberglass -reinforced plastic while the wing was made from magnesium 
to reduce its mass to a minimum. A three-phase induction motor was used 
to rotate the control surface in near sinusoidal oscillations about its 
midchord by means of the drive linkage shown in figure 2. The model was 
sting mounted in the wind tunnel and a strut was used to stiffen the body 
against the oscillating loads produced by the control surface. The loca- 
tion and dimensions of the strut are given in figure 1. 


(p = It - 2k radians 


(13) 


APPARATUS 


Wind Tunnel 


Model 
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There were four quantities obtained during this investigation: the 
amplitude of the oscillating lift induced on the wing in phase with the 
angular position of the control surface^ the amplitude of the oscillating 
lift induced on the wing in phase with the angular velocity of the control 
surface, and their respective centers of pressure. The measurement of 
unsteady forces such as those encountered in this investigation usually 
necessitates equipment in the nature of an oscillograph or an oscilloscope 
to record or view the dynamic signals. The use of such equipment is some- 
what inconvenient compared with that utilized in static tests in that, in 
many cases, the desired quantities are not indicated directly and addi- 
tional steps such as film reading and harmonic analysis must be added to 
the data-reduction process before the final result is obtained. It would 
be desirable, then, to use a technique that would eliminate these addi- 
tional steps and make it possible to read the peak magnitudes of the 
dynamic quantities directly. 

The technique and instrumentation employed in this investigation per- 
mitted the rapid, direct measurement of the peak magnitudes of the in-phase 
and out-of-phase components of the lift on the wing. The method consisted 
of multiplying a signal proportional to the induced lift on the wing by a 
sine or cosine signal in phase with the control-surface position and velo- 
city, respectively. These two products can be written as: 

Psin = sin u)t Lisin(o)t + cp) (l^) 

Pcos = cos wt Lisin(u)t + cp) (l 5 ) 


where 

Lisin(u)t + cp) . 
T 


Psin 

Pcos 

Equations ( lU) 


the oscillating lift induced on the wing 

phase angle between the oscillating lift and the angular 
position of the control surface 

product of the sine signal and the oscillating lift 

product of the cosine signal and the oscillating lift 


and (15) can be rewritten as: 



Li . / 

Psin = ~ [cos cp - cos( 2 u)t + 

cp)] 

( 16 ) 

Li . / 

Pcos = ~2 + 

Cp)] 

(17) 
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But the magnitude of the in-phase lift is defined as Licos 9 and that 
of the out-of-phase lift is defined as L^sin 9. Thus, it is evident 
that if the time-varying portions of the product signals are filtered out, 
the remaining parts are proportional to the components of the lift in phase 
with the angular position and angular velocity of the control surface. 

The technique described above can be realized with a system made up 
of three parts: a load measuring device capable of producing a signal pro- 
portional to the induced lift of the wing, a circuit producing the sine 
and cosine signals, and an indicating device that will respond only to 
steady signals. The force induced on the exposed portion of the wing was 
indicated by strain gages mounted on the wing supports. The arrangement 
of these strain-gage members is shown in figure 3* all, eight strain- 

gage bridges were used, two on each of the four wing supports. In order 
to compensate for any spanwise movement of the center of pressure, each 
Bridge on the left side of the body center line was connected in parallel 
with its counterpart on the right side. Thus there were four bridge pairs, 
two located at the wing leading edge and two located at the trailing edge. 
This arrangement of strain-gage circuits provided duplicate indication of 
the components of the induced lift acting at the leading and trailing edges 
of the wing, and the dual circuits were then used to find the in-phase and 
out-of-phase parts of these lifts. The s\jm of the in-phase lifts at the 
leading and trailing edges gave the total in-phase lift while the ratio 
of the lift at the trailing edge to the total lift gave the distance of 
the center of pressure from the leading edge in percent chord. These cal- 
culations were also carried out for the out-of -phase lift. 

The product of the force signal with a sine or cosine signal was 
obtained by using a sine or cosine voltage as power for the strain-gage 
circuits. The output of the gage circuit was then proportional to the 
stress in the gage multiplied by the applied voltage signal. The sine 
and cosine voltages were produced by a circuit which utilized a commercial 
induction resolver whose armature was driven by the motor driving the con- 
trol surface. A 20-kc sinusoidal carrier wave was used as an input signal 
to the resolver. This carrier wave was modulated sinusoidally by the 
resolver at the drive frequency to give two output signals with a 90° phase 
difference. These signals were each passed through a half-wave rectifier 
and filter to remove the carrier frequency and amplified to the level 
needed to drive the strain-gage circuit. 

The outputs of the strain-gage circuits were connected to light-beam 
galvanometers. The frequency response of these galvanometers can be 
described adequately by the response of a spring-mass -damper system. The 
differential equation of motion for the system is: 

mx + cx + kx = P sin ut (18) 

and the solution to this equation is 
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X 


^st 



sin(u)t + e) 


(19) 


where 


e 


tan“^ 



and 


Xst 


OJ 


o 


static deflection of the galvanometer to a steady signal of 
magnitude P . 


natural frequency of galvanometer. 



t damping ratio of galvanometer, — 

2 Vim 


Examination of equation (19) shows that if the ratio of the signal 
frequency w to the natural frequency Wq is sufficiently large, the 
response of the galvanometer to the signal is essentially zero. The 
response is also decreased by a high damping ratio. The galvanometers 
used during this investigation had sufficiently low natural frequencies 
and high damping that no response was observed to signals of 10 cycles 
per second or greater. As a result, the galvanometers acted as filters 
to remove the time -varying portions of the strain-gage signals and only 
the steady part of the signal was indicated. 


TESTS 


The investigation of the unsteady lift induced on the wing by an 
oscillating canard control surface was conducted over a Mach n-umber range 
of 0.6 to 0.8 and 1.4 to 1.9 at angles of attack of 0 , 5 , and 10 . The 
minimum supersonic Mach number was chosen such that the reflection from 
the tunnel walls of the shock waves from the control surface fell behind 
the wing trailing edge. A constant Reynolds number of 1. 6^X10 based on 
wing chord was held for the entire test. 

The control surface was driven at frequencies from 10 to 110 cycles 
per second covering reduced-frequency ranges of 0.06 to 0.70 at a Mach num- 
ber of 1.9 and 0.15 to 1.80 at a Mach number of 0.6. The maximum amplitude 
of the control surface oscillations was fixed by the drive linkage at 5°» 
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CORRECTIONS TO DATA 


The only corrections made to the data obtained in the subsonic speed 
range were to the Mach number and dynamic pressure of the flow. It was 
recognized that the phenomenon of wind-tunnel resonance and its effect 
on the air forces acting on an oscillating two-dimensional wing has been 
predicted on a theoretical basis and that experimental verification of 
this theory has been made (ref. 9)* However, the author knows of no pub- 
lished corrections for the effects of tunnel walls on the oscillatory air 
forces of a subsonic three-dimensional wing. Woolston and Runyan in ref- 
erence 10 give an equation for the resonant frequency of a tunnel having 
a rectangular cross section which has been used to calculate tunnel reso- 
nant frequencies for this investigation. These frequencies are shown in 
their appropriate places with the data and will be discussed more fully 
in a subsequent section. 

The Mach number and dynamic pressure of the subsonic flow were cor- 
rected for blockage effects by the method presented in reference 11. At 
a Mach number of 0.8, these corrections amounted to an increase of about 
1.5 percent in the Mach number and dynamic pressure over the value obtained 
from measurements made without a model in the tunnel. 

No corrections were made to either the flow conditions or the lift 
data at supersonic speeds. 

The possibility of errors due to air-stream fluctuations of a periodic 
nature was considered. These periodic fluctuations could be caused by 
tunnel turbulence, an oscillating wake shed from the nose of the model at 
angle of attack, or possibly by the wake from the supporting strut. In 
order to evaluate this effect, a series of tests were made on the model 
with the control surface removed. These tests showed that there were no 
periodic fluctuations in the air stream of measurable magnitude for the 
range of Mach mombers, frequencies, and angles of attack covered in this 
investigation. 


PRECISION OF DATA 


In the section "Technique and Instrumentation," the system used in 
this investigation was presented in its simplest form in order to show the 
essential features of its operation. As could be expected, actual practice 
proved the system to be more complex than was indicated from simple con- 
siderations. These complexities, as it turned out, determined the accuracy 
of the final results. In this section, the various factors that contribu- 
ted to a departure from the elementary system mentioned previously will be 
discussed in terms of their relationship to the precision of the data. 
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The usual starting point for a discussion of the precision of the 
data is concerned with the least readings of the quantities involved in 
the reduction of the data to coefficient form. For this investigation 
these readings were as follows: lift^ O.O5 pound; stagnation temperature^ 

2^ Fahrenheit; stagnation pressure^ 0.2 centimeter of mercury; period of 
oscillation (l/f)^ 0.0001 second. The total uncertainty in the coeffi- 
cients due to least readings was taken as the square root of the sum of 
the squares of the effects of these least readings. This calculation led 
to uncertainties of ±0.003 for ±0.12 for Opg^ ±0.005 center of 

pressure^ and ±0.5^ for cp.i 

In the earlier description of the technique^ the expressions for the 
terms in the product signals (eqs. (lU) and (15)) were written in an ide- 
alized form. In general^ however^ the strain-gage output cannot he 
expressed as a simple sine function but will also contain a steady term 
and harmonics of the fundamental frequency. These harmonics could arise 
from the distortion in the control-surface motion produced by the drive 
mechanism or from a nonlinear response of the downwash to the control 
deflection. The driving voltages could also contain harmonics introduced 
by the electronic circuitry producing them. Therefore^ the signals must 
be written as: 


00 

L - Lq +^LnSin(nu)t + cp^) 


(20) 


n=i 


^sin 


00 

= ^Bm sin(nia)t + 6m) 
n=i 


( 21 ) 


cos 


00 

cos(mwt + €ui) 


m=i 


where 

L total lift on the wing 
vsin sine voltage 
^cos cosine voltage 

^It should be noted that the uncertainty due to least reading becomes 
very large as the lift forces approach zero. The values given here are 
for average values of in-phase and out-of -phase lift. 
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Lq steady lift on the wing 

Ln amplitude of the nth harmonic of the oscillating lift 

q)j^ phase angle between the nth harmonic of the lift and the position 

of the control surface 

Bm amplitude of the mth harmonic of the sine voltage 

0JH phase angle between the mth harmonic of the sine voltage and the 
position of the control surface 

Cm amplitude of the mth harmonic of the cosine voltage 

£m phase angle between the mth harmonic of the cosine voltage and 

the velocity of the control surface 

The products of the lift and the sine or cosine signal are then 


00 ^ 

= BjiiSin(mwt + y Lf^BjnSin(nu)t + cpj^) sin(inu)t + 


sin 


m=i 


n=im=i 


00 00 TO 

= + 0^) + ^ |cos[(n - m)a)t + cpn “ " 


m=x 


n=i m=i 


cos[(n + m)u)t + cp^+ 0m] 


( 22 ) 


00 00 


Pcos = CmCOs(mu)t + e^) ^ LnCmsin(nwt + cpj^)cos(Tmot + e^) 

111=1 n=im=i 


00 00 00 

= ^Cjj^cos(mu)t + e^) j^sin[(n - m)u)t + cp^ - eml 


m=i 


n=i m=i 


sin[(n + m)wt + 



(23) 
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Examination of equations (22) and (23) shows that the only terms in 
the expressions that are not functions of time arise from the products of 
like harmonics in the lift signal and driving voltage. Thus the signals 
indicated by the galvanometers are: 

Psin = cos((pi - 0i) + cos(<p2 - 62) + • • • (24) 

Pcos = sin(cpi - ei) + sin(cp2 - €2) + • . • (25) 

The desired quantities are, of course, the in-phase and out-of -phase com- 
ponents of the fundamental amplitude of the lift, that is, Licos (pi and 
Lj^sin Tx • 

It can be seen then that the higher harmonics can produce spurious 
signals which would lead to erroneous results. 

Since no measurements were made of the harmonic content of the lift 
signal, it was not possible to determine quantitatively what errors were 
introduced from this source. However, harmonic analyses of the output of 
the resolver circuits were made during calibration and it was found that 
there was 7 percent of the fundamental amplitude present as second harmonic 
and 1 percent as third harmonic in both the sine and cosine voltages for 
the frequency range of the investigation. All higher harmonics were negli- 
gible. Thus if Bi = 1, then B2 = O.O7, B3 = 0.01 and B4, B5, . . . = 0. 
Hence, if the amplitude of the second harmonic equals that of the funda- 
mental in the lift signal, the error in the measurement of the fundamental 
amplitude would be 7 percent. It can be seen therefore that the low har- 
monic content of the sine and cosine voltages reduces any error due to 
harmonic content of the lift force. 

The misalinement of the resolver was another factor that could lead 
to error in the final result. This misalinement can be thought of in terms 
of a phase angle between the sine and cosine voltages and the position and 
velocity of the control surface, respectively. The resolver would be per- 
fectly alined if these phase angles were zero. The general form of the 
products written above incorporated these phase angles. Inasmuch as it 
has been shown that the effects of the second harmonics, or higher, are 
small, only the alinement of the fundamental amplitude of the resolver 
signals was considered. Calibrations were made during the course of 
testing which showed that the average values of 0i and ex were about 
0.5° each. It can be shown that a resolver misalinement of 0.5° will 
introduce an error of not more than 0.5° in the calculation of the phase 
angle (Pi . 

The remaining factor to be discussed with regard to the precision of 
the data is the response of the mechanical system to oscillating loads. 
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It was recognized that the output of the strain gages indicating the lift 
on the wing would he influenced by the response of the various parts of 
the mechanical system to the oscillating loads. In order to obtain a 
quantitative measure of this influence, a series of calibrations were made 
on the model in place in the wind tunnel by applying an oscillating force 
of known magnitude and measuring the amplitude and phase response of the 
system indicated by the output of the strain gages as the frequency of the 
input force was varied. These calibrations showed that above a frequency 
of 110 cycles per second, the natural frequency of the wing panels in 
bending (l 40 cps) introduced large errors in both amplitude and phase 
angle. In addition, it was found that the body natural frequency in bend- 
ing occurred at approximately 85 cycles per second. However, the influ- 
ence of this resonant peak on the response of the system was restricted to 
a narrow frequency band. Thus, by testing over a frequency range up to 
110 cps and omitting a band near 85 cps, large errors arising from reso- 
nance of parts of the mechanical system were eliminated. With these 
regions eliminated, the results of the calibrations still showed small 
deviations from the ideal response. Since this scatter was greater than 
any known inaccuracies in the instrumentation, it was felt that it should 
be attributed to the response of the mechanical system and therefore had 
to be taken into account in the precision of the results. The calibra- 
tions showed that the root-mean- square deviation of the indicated force 
from the true force was ± 6.4 percent of the true force, the rms deviation 
of center of pressure was ±0.02 wing chord, and the rms deviation of the 
phase angle from zero was ±1.7^* 

It can be seen that the scatter in the results of the frequency- 
response calibrations for the mechanical system is the primary factor in 
the determination of the precision of the final result. The calculation 
of the uncertainty in the results including the effects of least reading, 
resolver misalinement, and mechanical-response errors led to values ± 0.024 
for Cl , ± 0.24 for Cl^, ±0.02 for x/c, and ±1.9^ for cp. 


RESULTS AKD DISCUSSION 


The results of the experimental investigation of the unsteady lift 
induced on a wing in the downwash field of an oscillating canard control 
surface are shown in figures 4 through 8 as a function of reduced frequency 
for several Mach numbers and angles of attack. The nondimens ional lift 
derivatives in phase and out of phase with control-surface position are 
presented in figures 4 and 5, while their respective centers of pressure 
are shown in figures 6 and 7. The variation of phase angle with reduced 
frequency is shown in figure 8 for the data obtained at a = 0 ^. It should 
be noted that, since the experimental values of the lift derivatives were 
obtained by dividing the lift by ^o^ assumption is implicit that 

lift is a linear function of 5 and 5 . 
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Also shown on these figures are the corresponding theoretical values 
computed using equations (lO)^ (H)^ 3*nd (l3)- values of static lift, 

Cl^*, used in these calculations were computed hy the method of refer- 
ence 1. Thus, while the theoretical development of equations (lO), (ll)^ 
and (13) was based on a body less model, the calculated values of 
were based on wing -body theory. In other words, it has been assumed that 
the body affects the initial values of the lift derivatives (k-^O) but 
does not affect their dependency upon frequency given in equations (lO) 

and (11). 

In-phase lift derivative . - Figure ^ is a plot of the in-phase lift 
derivative versus reduced frequency for the Mach numbers and angles 

of attack covered. Examination of this figure shows that the theory of 
reference 1 provides a reliable guide to the estimation of at low 

frequencies for an angle of attack 0^. In addition, the theoretical trends 
of Cpg with frequency calculated using the simplified indicial response 

were, in general, borne out by the experimental data for this angle of 
attack. At the higher angles of attack, however, there were some depar- 
tures from the theory. For supersonic Mach numbers, the primary differ- 
ences were in magnitude. The trend of the data with frequency remained 
the same as the angle of attack increased. At subsonic speeds, large 
variations in both magnitude and trend with frequency appeared. This was 
most noticeable at an angle of attack of 10*^. While there is not suf- 
ficient evidence to explain this result conclusively, partial separation 
of the flow over the control surface may be one of the causes. Static 
tests of a similar wing (ref. 12) show that at subsonic speeds, the flow 
begins to separate at an angle of attack of about 10^ and that the lift no 
longer increases with increasing angle of attack above l6^. Since the 
angle of attack of the control surface varied from 5^ "to 15^ during each 
cycle when the body angle of attack was at 10^, the character of the flow 
was probably changing radically through each cycle, thereby affecting the 
lift Induced on the wing. 

The possibility of an effect of wind-tunnel resonance on the data 
obtained at subsonic speeds must be considered. The wind-tunnel resonant 
frequencies, k^, are shown in the figures in terms of reduced frequency 
for each of the subsonic Mach numbers. It was shown in reference 9 that, 
for the two-dimensional case, the effect of tunnel resonance is to decrease 
the lift markedly near the resonant frequency. Examination of the data 
obtained in this investigation indicates that if a wind-tunnel resonance 
phenomenon exists for this particular combination of model and tunnel, it 
does not have the simple effect indicated in reference 9- 

A comparison of the results calculated by the simplified indicial 
response and those calculated by the more exact indicial response is shown 
in figure 4(b) for a Mach number of 1.7- It is apparent that the differ- 
ences are quite small over a large portion of the frequency range, thus cor- 
roborating the statement concerning this fact made in section "Theory.” 


CONFIDENTIAL 


18 


• • • • 
• • • 




NACA RM A55P01 


It may be well to mention here the point brought out in the intro- 
duction concerning artificial damping for canard-type missiles. Since 
damping may be added by motion of the control surface proportional to the 
missile pitching velocity, the in-phase lift shown in figure 4 contributes 
to the damping of the missile. The results show that if the pitching 
frequency is sufficiently high, there would be a serious reduction and 
possibly a reversal of the damping moment supplied by the rear surface. 
However, the probability that the pitching frequency will reach these 
values appears small. 

Out-of -phase lift derivative . - The lift derivative out of phase with 
control deflection, is shown plotted versus reduced frequency in fig- 

ure 5 . Here the predictions of the theory are confirmed by the experiment 
for a = 0^ at subsonic speeds. However, as was the case for the in-phase 
derivative, the agreement between experiment and theory became poorer as 
the angle of attack increased, with the data obtained at a = 10° showing 
a marked difference from that obtained at the lower angles. 

In the supersonic speed range, the theoretical values of at 

low frequencies show somewhat better agreement with experiment at a = 0° 
than ab the higher angles of attack. However, the decrease in with 

increase in the frequency predicted by the simplified indicial response 
is not reflected in the experimental values. Examination of figure 5('t>) 
shows that the more exact indicial-response result is similar to that for 
the simplified indicial response. Calculations of the values of 

have shown that the frequency response of this derivative is highly sen- 
sitive to changes in the details of the indicial curve. Hence, better 
agreement could be achieved with a more refined indicial function. 

Center of pressure of in-phase lift .- The experimental values of the 
center of pressure of the in-phase lift are shown in figure 6 along with 
the theoretical values for a = 0°. The theoretical values were calculated 
by the method of reference 1. Since this method is based on linear theory, 
no variations of center of pressure with angle of attack are predicted. 

The theory of reference 1 predicts the location of the center of 
pressure very well for all Mach numbers at low frequencies and low angles 
of attack. In addition, the lack of variation with frequency predicted 
by the simplified indicial approach is borne out by the experimental data 
over most of the frequency range. The large deviations from the theoreti- 
cal values were obtained either for conditions where the in-phase lift was 
near zero, thus affecting the accuracy of measurement, or at a = 10^ for 
subsonic Mach numbers where it is known that the lift forces themselves 
do not follow the predicted results. 

Center of pressure of out-of -phase lift . - The experimental data for 
center of pressure of the out-of-phase lift, presented in figure 7^ 


CONFIDENTIAL 


NACA EM A55F01 •*: \V :*‘Ci)NF^DpTj:*A]5|* 

••• • • ••• • • 


19 


confims, in general, both the magnitude and trend with frequency pre- 
dicted by the theory (excluding the data at a = 10^) . At the lower 
frequencies, the data obtained in the subsonic speed range show values 
somewhat forward of those predicted by theory, but at supersonic speeds 
agreement is good throughout the frequency range. 

Phase angle . - Perhaps the most direct application of the lag-in- 
downwash concept is to the calculation of the phase angle of the lift 
induced on the wing. Here, the time necessary for a disturbance to travel 
from the foirward surface to the rear surface - the lag - can be converted 
directly to a phase angle for a given frequency. Figure 8 shows a com- 
parison between this theoi*y and the phase angles obtained at all Mach 
numbers for an angle of attack of 0°. It can be seen that this simple 
theory predicts the experimental variation of phase angle with frequency 
adequately for all Mach numbers for values of k up to about 1.0. Beyond 
that value the experimental results are somewhat higher than theory. While 
only the phase angles for a = 0° are presented in this figure, the phase 
angles for the other angles of attack, with the exception of those obtained 
at a = 10° and subsonic speeds, showed the same variation with reduced 
frequency. 


CONCLUSIONS 


An experimental investigation of the unsteady lift induced on a wing 
by the downwash field of an oscillating canard control surface led to the 
following conclusions: 

1. Existing theories provided a reliable guide to the estimation of 
the magnitudes of the in-phase and out-of-phase lift derivatives and their 
respective centers of pressure at low values of reduced frequency and low 
angles of attack. 

2. The trends of the data with frequency calculated with the use of 
a simple indicial lift response for the wing were, in general, confirmed 
by the experimental results at low angles of attack. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., June 1, 1955 
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Figure 1.- Sketch of model. 
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Figure 2 .- Photograph of drive system. 
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Figure 3*“ View of wing and strain-gage members. 
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(a) Subsonic Mach numbers. 

Figure 4.- Variation of in-phase lift derivative with reduced frequency. 
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(b) Supersonic Mach numbers. 

Figure 4 .- Concluded. 
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(a) Subsonic Mach niuabers. 

Figure 5.- Variation of out-of-phase lift derivative with reduced frequency. 
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(b) Supersonic Mach numbers. 
Figure Concluded. 
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(a) S-alDSOnic Mach minibers. 

Figure 6.- Variation of center of pressure of in-phase lift with reduced 

frequency. 
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(b) Supersonic Mach numbers. 
Figure 6.- Concluded. 
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(a) Subsonic Mach numbers. 

Figure 7.- Variation of center of pressure of out-of-phase lift with 

reduced frequency. 
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(b) Supersonic Mach numbers. 

Figure 7«- Concluded. 
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Figure 8 .- Variation of phase angle with reduced frequency; 


a = 0^ 
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